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THE PROBLEM OF INTERNAL AND EDGE CRACKS 
IN AN ORTHOTROPIC STRIP* 

by 

F. Delale and F. Erdogan 
Lehigh University 


ABSTRACT 

The plane elastostatic problem of internal and edge cracks in an 
infinite orthotropic strip is considered. The problems for the ma- 
terial types I and II are formulated in terms of singular integral 
equations. For the symmetric case the stress intensity factors are 
calculated and are compared with the isotropic results. The results 
show that because of the dependence of the Fredholm kernels on the 
elastic constants, unlike the crack problem for an infinite plane, 
in the strip the stress intensity factors are dependent on the elastic 
constants and are generally different than the corresponding isotropic 
results. 


*This work was supported by NASA-Langley under the Grant NGR-39-007-011 
and by NSF under the Grant ENG 73-045053 A01 . 
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1 . INTRODUCTION 


In plane elastostatic problems for an infinite orthotropic medium 
containing a line crack [1-3] or a series of collinear cracks [4] it 
was shown that the stress intensity factor is identical to that found 
for isotropic materials. However, if the geometry of the medium is 
bounded, it is expected that in orthotropic solids the material con- 
stants would influence the stress intensity factors. A bounded speci- 
men geometry which is sufficiently simple for the purpose of analysis 
and at the same time is of sufficient practical interest is that of a 
long strip containing internal or edge cracks. The main objective of 
the present paper is by considering this problem to give some idea 
about the degree of influence of the material orthotropy on the stress 
intensity factors. The equivalent isotropic case is one of the more 
widely studied crack problems in technical literature {see, for ex- 
ample, [5-13]). 

2. FORMULATION OF THE PROBLEM 


Consider the plane problem for the orthotropic strip shown in 
Figure 1. Referring to, for example [14], the equilibrium equations 
for an orthotropic plane may be expressed as 



9 z u + a 
Sy 2 " p 3 9x9y 


= 0 



+ 3 


9 Z V 
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9 z i 


2 3y^ 3 9x9y 


= 0 


( 1 ] 


where u,v are the x,y-components of the displacement vector and 
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for generalized plane stress, and 
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for plane strain. Here, E^, v^, are the engineering elastic con- 
stants, (i , 0 ) « (1,2,3), the indexes (1,2,3) corresponding to the direc- 
tions (x,y,z), and the matrix (b..) is given by 

* J 


(b i j) = B = A” 1 , A=(a ij( ) , (i,j) = (1,2,3) 
a ii = 1/E ii 5 a ij = " V ij /E ii = a ji * (i?4j) 


(4) 


The solution of the problem shown in Figure 1 may be obtained by using 
the standard superposition technique. Thus, for the purpose of evaluat 
ing the stress intensity factors and obtaining information relevant to 
the fracture of the solid, it is sufficient to consider the problem in 
which statically self-equilibrating crack surface tractions are the 
only external loads. 


To solve the differential equations (1) let 


u(x,y) «* ~ [f-j (a.x)cosay + g 1 (a,y)sinax]da , 
v(x,y) = j J [f 2 (a,x)sinay + g 2 (a,y)cosax]dct . 
Substituting from (5) into (1) the functions f. and g.. , (i=l,2) are 


,00 


(5) 


obtained as follows: 


f-j (a>x) = l A..(a)e J , f 2 (a,x) = l CjAj(a)e J , 

(6 

“ s.ay/05 “ s-ay/Ps 

9-j (a, y) = J Bj(a)e J , g 2 (a,y) = J d,.Bj(a)e J , 

where are the roots of the following characteristic equation 


s * + 8 4 s 2 +p 2 = 0 , s 3 = -s-j , s 4 = -s 2 , 
and the constants 3 C , c., and d., (j=l,..,4) are given by 

4 J J 

P 4 = (3| - &1&2 " ^/^i » B 5 = B 2^1 ' 

C-J = ~c 3 ~ (1 ~ » c 2 ~ -c 4 ~ ( 1 “ 3i s 2^3 s 2 * 

d] = -d 3 - (s| - , d 2 = “d^(s 2 - ^|)/^3 S 2^5 


Assuming that x and y are axes of symmetry for loading as well as 
geometry, the unknown functions A. (a) and B.(a), (j=l,..,4) are de- 

J J 

termined from the following conditions: 


u(x,y)+0 , v(x,y)-K) for y*» , 


tf xx (h,y) = 0 , o xy (h,y) = 0 , 0<y<°° , 


u(0,y) = 0 , a xy (°»y) = 0 > 0<y<~ , 
° xy { x »0) = 0 > 0<x<h , 


°yy( x,+ °) = " p M 5 a< M <b * 


v(x,0) = 0 , 0< | x | <a , b< j x|<h , j 



where the crack surface traction p(x) is a known function. The seven 
homogeneous conditions (9-12) may be used to eliminate seven of the un 
knowns and the mixed boundary conditions (13) would give a system of 
dual integral equations to determine the eighth. In this paper, de- 
fining a new unknown function 

<J>(x) = —• v(x,0) , a<jx|<b , (14 

the problem is reduced to a singular integral equation by using the 
first equation of (13). From the second equation of (13) it is seen 
that 

<|>(x,0) = 0 , 0<|x|<a , b< | x l <h , 

rb 

L <}>(x)dx - 0 (IE 

Examining the roots of (7) it may be observed that 

(i) for 3 4 <0, 3 6 = 3|-43g>0 there are four real roots, s-j , s 2 
Sg—Sp and s 4 = -s 2 (s-j>0, s 2 >0); in this case the corresponding 
material will be classified as type I; 

(ii) for 3 6 <0 the roots are complex, s^ = w-j + iw 2 = -s 3 » 
s 2 = ^-iw 2 =-s 4 (^>0, w 2 >0) and the related material will be classi- 
fied as type II; and 

(iii) for 3 4 >0, 3 g >0the roots would be pure imaginary 
S 1 ~ ^ w 3 “ -s 3 * s 2 = ^ a) 4 = " s 4‘ 

In practice generally 3 4 is negative. Therefore the problems of 
interest are those relating to materials type I and II only. This 



classification seems to be necessary in order to pursue the formulation 
of the problem beyond equations (6) without introducing unnecessarily 
complicated complex algebra. Also, in this paper only the case of 
generalized plane stress will be considered. For plane strain it is 
sufficient to replace the quantities E x /0-^ xy v yx ) » E y /(l-v xy Vy X ) , and 
v yx E x^ 1_v xy v yx^ by b ll » b 22’ and b 12’ res P ective1 y (see equation 4). 

Because of symmetry considering only one quarter of the medium 
shown in Figure 1, and using the standard stress-displacement relations 
for plane stress, after somewhat lengthy but routine analysis, for ma- 
terial type I (i.e., for real s-j and s 2 ) the problem may be reduced to 
the following integral equation: 


f C^ + t?x +k i{ x ’t) " k-j (x s -t)3<f-(t)dt 

^ cl 

41-v v ) 

= 2Eyr^“ > a<x<b > 

under the additional condition (15). Here the kernel is given by 

-.(h-t)aB s /s 1 


(16) 


¥ Xit) = s^|o [K i (x >“ )e 

-(h-t)ap 5 /s 2 


m 14 
+ K 2 (x,a)e 


]da . 


(17) 


The expressions for K-j, Kg, and m-j 4 are given in Appendix A. 

For material type II the roots of the characteristic equation (7) 
are complex. Defining 


V V ia) 2 = ' s 3 


s 2 = a>i 


•™ 2 = ~s 4 


(18) 
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and, assuming that w^>0, in this case the integral equation becomes 
f b [^+ tlx+ k 2 (x,t) - k 2 (x,-t)]*(t)dt 

i 3l 

1 T ( 1 " V vAv) 

= “ PE r y ' y , X PM * a<x<b » 09) 

ZE y r 14 

again subject to condition (15). The kernel k 2 is given by 

(“ -w,a(h-t) 

k 2 (x,t) = I K 3 (x,t,a)e da , (20) 

where the function K 3 and the related constants r^ are defined in the 
Appendix B. 


One may note that in the special case of single internal crack, 
(i.e., for a=0, b<h) the integral equations (16) and (19) may be 
written as 


/■b i 

| [^17+ k i (x,t)](j>(t)dt “ - jjj7 p(x) 
-b 1 


-b<x<b 


(i=l,2) , M-j s i , M 2 = l-v Y v 

xy yx xy yx 


( 21 ) 


where i=l and 2 correspond to materials type I and II, respectively. 
In this case the single-valuedness condition (15) becomes 


fb 

<|>(x)dx = 0 . 

ib 


( 22 ) 


3. STRESS INTENSITY FACTORS 


The standard definition of the stress intensity factors at the 
crack tips a and b is 

k(a) = 11m /2U-x; a (x, 0) , 

x+a yy 

k(b) = Tim ojx.0) . (23) 

x*>-b yy 

To determine the asymptotic behavior of the cleavage stress around 
the crack tips, first it may be observed that the index of the singular 
integral equations (16) and (19) is +1. Therefore, the solution is of 
the following form [15]: 

<j)(t) = f(t)[(t-a) (b-t)] . (24) 

Next, it should be pointed out that the left hand side of (16) and (19) 
gives a (x , 0 ) for x outside the interval ( a , b) as well as within. 

J J 

Thus, making a change in variable 

x = + jt = ^5. r + b+a (25) 

for example, (16) may be expressed as 

i|\i + ^ Sjr)] , (r )d r =^ , b±|< s <21l^za (26) 

where 

q(s) = c (x,0) , i|>(r) = <l>(t) = F(r)(l-r 2 ) 


(27) 




and k(s,r) and F(r) are bounded functions. The objective is then to 
determine the asymptotic behavior of q(s) around s-+l, (|s|>l) In terms 
of the unknown function F(r) which is obtained by solving the Integral 
equation (26) in -l<s<l where q(s)=-p(x) is known. To do this let us 
assume that the bounded function F(r) can be represented by the follow- 
ing infinite series in Tchebyshev polynomials: 

F(r) = 2 A n T n (r) * (28) 

0 

Substituting from (28) into (26) one obtains 


lA„[e„(s) + H n (s)] 


g(s) 

Mi 



T n ( r ) dr 

(r-sJ/Tr 7 


H n (s) = 


1 

IT 


k(s»r)T n (r) 

1 1 “TPF 7 


dr . 


Here, H n (s ) is a bounded function. To determine G n one may use 
expression 


(29) 

(30) 


the 


1 f 1 T n (r)dr = (z-v?CD n 
B !i (r-z )/prr “ 


(31) 


where z is the complex variable in the plane cut along (-1,1). Ob- 
serving that on the real line (z 2 -l is an odd function, from (30) 
and (31) it follows that 

G ( s ) = - Is-sgnlsbgZP!!. . (32) 

n sgn(s)v^T 

As s-41 (32) yields 
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s*+l» (s>l): G n (s) = - + R 1 (s) , 

s+-1. (s<-l): G n (s) = + R 2 (s) , (33) 

where the functions R-j and Rg are bounded at s=+l. 

Now, observing that H n (+1) = finite, T n (l) = l, T n (-1) = (-1 )\ from 
(29) and (33) the asymptotic behavior of q(s) is found to be 

s++l, (s>l): - -^j=+ Rj(s) , 

s+-l. (s<-D: ^j' + V s > ' (34) 

where again the functions R 3 and R^ are bounded at s=+l. Going back 
to the original quantities by using (25) and (27), (34) becomes 

x*b, tx>b): = . .am + R6 ( X) , 

XH-a, (x<a): + r 6 ( x ) , (35) 

where the functions R g and R g are also bounded at x=b and x=a. Thus, 
from (23) and (35) the stress intensity factors are found to be 

k(b) = - I^FOMT , k(a) = M^-l)/* . (36) 

In the case of fully imbedded cracks the integral equations (16), 
(19) or (21) can always be reduced to the normalized form (26) and can 
be solved by using the technique described in [16]. 



4. EDGE CRACKS 

In equations (16) and (19) the kernels k-|(x,t)-k-j(x,~t) and 
k 2 (x,t)-k 2 (x,~t) are bounded provided b<h (see Figure 1). For b=h, 
that is in the case of edge cracks, the integral equations are still 
valid but these kernels do not remain bounded as x and t go to the 
end point b=h and, consequently, the singular behavior of the solu- 
tion at x=b=h may no longer be described by (24). Expressing the 
kernels in (16) and (19) as 

k.j(x,t) s k. f (x,t) + k. s (x,t) , (1=1 >2) 

where k^ is bounded in the closed interval [a,h], the unbounded parts 
k is , (1-1 ,2) may be obtained from (17) and (20) by examining the asymp 
totic behavior of the integrals for large values of a. Thus, after 


some routine analysis we find 








( 38 ) 



rggOJjUh-t-x) 

+ uiflSh-t-xr + ^Uli'-t-xP' 


where the constants m. and are given in the appendices, 
example, the integral equation (16) may be expressed as 

| Ct^ + Fx +k ls (x>t) + k if( x>t ) " k i(x/t)]<|»(t)dt 
= - J p(x) , a<x<h . 


Thus, for 


(39) 


In (39) for the purpose of asymptotic analysis transferring the 
terms involving the bounded kernels to the righthand side one may 
write 


I' EtV k ls (x.t)]*(t)dt = P^x) , a<x<h 

where P-j is a bounded function in [a,h]. Letting now 
f,(t) 

<J>(t) = £ y , 0<Re(a,p)<l , 

(h-t) a (t-a)^ 


(40) 


(41) 


where f-j is H-continuous in [a,h], and following the procedure outlined 
in [15], the characteristic equations for a and B are found to be 


cotnB = 0 , B = 1/2 , (42) 

-costo + - ~ - [m^tBs/sf^fSi/Bs) + m 17 (si/B 0 ) 

14 * 15 

+ n> 18 (s 2 /B 5 ) + m l9 (B 5 /s^) a (s 2 /B 5 )] = 0 . (43) 


- 12 - 


Similarly, for material type II described by (19), assuming the 
solution again as given by (41), the characteristic equations become 

cotTrp = 0 , 0 = 1/2 (44) 


cosira + 


_J 

r 14 r 19 


r 29 


+ 


a Z r ZS~^ r Z 7 


cos(2atan 



^1 r 25' fa) 2 r 27 
on + wi 


-1 w 2 

sin(2atan ^-) 


(45) 


At the imbedded crack tip x-~a it is seen that the singularity has the 
expected 1/2 power. On the other hand, as in the isotropic case, (43) 
and (45) have no root for which 0<Re(a)<l, meaning that at x=h there is 
no power singularity. One may also proceed and investigate the possi- 
bility of a logarithmic singularity for the solution. Thus, letting 
a=0 in (41) and defining the sectionally holomorphic function 



we find [15] 


F,(z) 


simrg 


i , 

(z-a) e ir(h-a) 8 


log(z-h) + P 2 (z) , 



1 f b (t(t)dt 

tt J a t-(2h-xT 


f-j (a)cotirg 
<x-a) 8 

F 1 (2h-x) 


+ 


^(h) 

ir(h-a)^ 


log(h-x) + P 3 (x) , 


(47) 


9 • • * 


where Pg and Pg are bounded at x=h and have at most a singularity of 
lower order than 0 at x=a. Substituting from (47) into (40), multi- 
plying through by (x-a)^ and letting x+a, it is found that cotir0=O, 
giving again 3=1/2. On the other hand, in the neighborhood of the 
end point x=h one obtains 

O -m^j£ (5,111,6+5,111, 7 +s 2 «. 18 +s 2 m 19 )]log(h-x) + P 4 (x).P 1 (x) (48) 

where P^ contains all the bounded terms around x=h on the left hand 
side of (40). Similarly, for the material type II one finds 

^ + '(w^)r u r lg ( V25“V26" w l r 27' w l r 28^ 1o ^ h “ x ^ + ^4 (x) =Q l (x) 

(49) 

It turns out that, as in the case of isotropic materials [10], the co- 
efficient of the logarithmic term in (48) and (49) is identically zero, 
meaning that the solution may not have logarithmic singularity at x=h. 
In the edge crack problem the integral equation (39) and the similar 
equation for the material type II are solved by defining 

$(x) = &A = JEM (50) 

/s+T 

and by using the numerical technique described in [10]. In this case 
the stress intensity factor at x=a becomes 


k(a) = M,F(-1)/2I , % - (h-a)/2 . 


(51) 


5. RESULTS AND DISCUSSION 


As an example the following two orthotropic materials will be 
considered: 

Type I: F. = 8x10^ psi (55.16 x 10 9 N/m 2 ) , 

A 

E y = 24.75xl0 6 psi (170.65 xlO 9 N/m 2 ) , 

G yv/ = 0.7 xlO 6 psi (4.83xl0 9 N/m 2 ) , 

Ajr 

V xy = 0.036 , 

Type II: E ='3.1 x 10 6 psi (21 . 37 x 10 s N/m 2 ) , 

A 

E y = 9.7 xlO 6 psi (66.88 x 10 s N/m 2 ) * 

& = 2.6 x10 s psi (17.93x10 s N/m 2 ) , 

Ay 

\y-0.Z . 

Tables 1-3 show some of the calculated results for the stress inten- 
sity factors, In all the calculations it was assumed that the crack 
surface traction was constant, l.e. > 

ffyyU.O) = -P(x) * -P 0 (52 

which corresponds to uniform tension of the strip away from the crack 
region. Table 1 shows the results for an internal crack of length 2b 
(see Figure 1, a=0) which was found by solving (21). The stress in- 
tensity factors used in the tables are defined by (23) and are calcu- 
lated by using (36) for internal and (51) for edge cracks. The stress 
intensity factors for symmetrically located two collinear cracks 
(Figure 1) are given by Table 2. Table 3 gives the results for sym- 
metric edge cracks. 




The tables also contain the stress intensity factors for the 
isotropic strip which are included for comparison. A close examina- 
tion of the integral equations (16), (19), or (21) would indicate that 
in orthotropic materials since the Fredholm kernel k-j or k 2 is heavily 
dependent on the material constants, the solution must also depend on 
the constants. On the other hand, in isotropic materials even though 
the structure of the integral equation is identical to that of (16) or 
(19) (see, for example, [10]), the kernel of the integral equation is 
independent of the elastic constants and the constants appear in the 
equation as a multiplying factor (in the form of (1 +k)/4p) only. The 
stress intensity factors given in the tables indicate that the results 
for the orthotropic strip are indeed different than the isotropic re- 
sults. The tables also show that for approximately the same modulus 
ratio Ey/E x (in the example approximately 3/1), depending on the re- 
maining constants, the materials may not only be of different type 
(I or II), also the stress intensity factors may be greater (in this 
case, in material type II) or smaller (in material type I) than the 
isotropic values. In orthotropic materials there are three independent 
material parameters, namely, G /E u , E /E. and v . Therefore, it 
does not seem to be feasible to make a systematic study and demonstrate 
the effect of the material orthotropy on the stress intensity factors. 
However, it appears that there exists a difference between isotropic 
and orthotropic results and in highly orthotropic materials it may 
be significant. 

In solving this problem, the numerical analysis produced a some- 
what unexpected result. First, it should be pointed out that the 
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results given in the tables are accurate to roughly three significant 
digits, the remaining digits may not be accurate. On the other hand, 
after rotating the material 90 degrees (i.e., taking the strip parallel 
to the less stiff axis and the crack along the stiffer axis) and fully 
expecting to obtain a different set of results, the print out for the 
stress intensity factors came out to be identical - in all eight digits ■ 
to the original values obtained for the strip which was parallel to the 
stiff axis. Furthermore, the ratio of the function F(r) defined by (27) 
at all points in -l<r<l for the 0 and 90 degree orientations was found 
to be constant, indicating that the crack surface displacements for the 
two cases are related by (see (14), (27), and (36)) 

v 0 (x,0)M o - v 90 (x,0)M 90 (53) 

where the constant M is defined by (21). This simply shows that the 
kernels k^ and kg which appear in the integral equations (16), (19), 
and (21) remain invariant under a 90 degree rotation for a given ortho- 
tropic strip. 

Table 1. The stress intensity factor k(b)/p 0 v^B' for an internal crack 
of length 2b in isotropic and orthotropic strips. 


U '/ L. 


Orthotropic 

D/n 

i soiropic 

Type I 

Type II 

>0 

+1.0 

+1.0 

+1.0 

0.1 

1.0060 

1.0044 

1.0064 

0.2 

1.0246 

1.0182 

1.0261 

0.3 

1.0578 

1.0428 

1.0611 

0.4 

1.1094 

1.0811 

1.1155 

0.5 

1.1869 

1.1387 

1.1966 

0.6 

1.3033 

1.2264 

1.3183 

0.7 

1.4888 

1.3674 

1.5099 

0.8 

1.8160 

1.6241 

1.8471 

0.9 

2.5809 

2.2487 

2.6278 



Table 2. The stress intensity factors k(a) and k(b) for symmetric 
collinear internal cracks in a strip. 


k(a)/p 0 /£ 


k(b)/p Q /£ 


0.4 

0.5 

0.6 

0.8 

0.9 

+1 

0.9 

0.95 

0.98 

*>1 


Isotr. 


{-HO) 


176 
109 
097 
127 

+1.600 
678 


Type I Type II Isotr. Type I Type II 


60 

95 


1.081 

1.104 

->1.531 

1.595 

1.160 

1.226 

+1.600 


+1.613 

1.689 

1.200 

1.275 

+1.661 


1.094 

1.122 

1.221 

M 

1.694 

1.445 

1.875 

(+») 


i . 1 1 
1.096 
1.127 
1.231 

(•**) 

1.705 

1.461 

1.883 

M 


Table 3. The stress intensity factor k(a)po/£ for the case of 
symmetric edge cracks, 2,= (h-a)/2, 


a/h 

k(a)/p Q /& 

Isotropic 

Type I 

Type II 

0.1 

2.980 

2.978 

2.982 

0.2 

2.218 

2.208 

2.220 

0.3 

1.907 

1.887 

1.912 

0.4 

1.742 

1.710 

1.750 

0.5 

1.640 

1.600 

1.661 

0.6 

1.600 

1.531 

1.613 

0.7 

1.574 

1.486 

1.590 

0.8 

1.567 

1.462 

1.587 

0.9 

1.576 

1.458 

1.593 

0.98 

1.582 

1.467 

1.598 





















REFERENCES 


1. D, D. Ang and M. L. Williams, "Combined Stresses in an Orthotropic 
Plate Having a Finite Crack", J. A ppl. Mee/t., vol . 28, Trans. ASME, 
pp. 372-378, 1961. 

2. G. N. Savin, Stress Concentration Around Holes, Pergammon Press, 

New York, 196T: 

3. G. C. Sih, P. C. Paris, and G. R. Irwin, "On Cracks in Rectilinearly 
Anisotropic Bodies", Ini. J. FAaciuAe MzcJianias, vol. 1, pp. 189- 
203, 1965. 

4. S, Krenk, "Stress Distribution in an Infinite Anisotropic Plate 
with Coll inear Cracks", Ini. J. Solid* S&matu/Le* , vol. 11, 

pp. 449-460, 1975. 

5. W. T. Koiter, "On the Flexural Rigidity of a Beam Weakened by 
Transverse Saw Cuts", Paoc,. Royal Nei/i. Acad, oh Scitnce*, B59, 
pp. 354-374, 1956. 

6. W. T. Koiter, "Rectangular Tensile Sheet with Symmetric Edge Cracks", 
J. A ppl. Mec k. f vol. 32, Trans. ASME, p. 237-242, 1965. 

7. H. F. Bueckner, "Weight Functions for the Notched Bar", ZAMM, 
vol. 51, pp. 97-109, 1971. 

8. H. F. Bueckner, "A Novel Principle for Computation of Stress In- 
tensity Factors", ZAMM, vol. 50, pp. 529-545, 1970. 

9. 0. L. Bowie, "Solutions of Plane Crack Problems by Mapping 
Technique", Methods of Analysis and Solutions of Crack Problems , 

G. C. Sih, ed., Noordhoff International Publishing, Leyden, 1973. 

10. G. D. Gupta and F. Erdogan, "The Problem of Edge cracks in an 
Infinite Strip", J. A ppl. M &ah., Trans. ASME, vol. 41, pp. 1001- 
1006,1974. 

11. I. N. Sneddon and R. P. Srivastav, "The Stress Field in the 
Vicinity of a Griffith Crack in a Strip of Finite Width", Ini. J. 
Engng. Sal., vol. 9, pp. 479-488, 1971. 

i 

12. M. Isida, "Laurant Series Expansion for Internal Crack Problems", 
in Methods of Analysis and Solutions of Crack Problems , G. C. Sih, 
ed., Noordhoff International Publishing, Leyden, 1973. 

13. S. Krenk, "On the Elastic Strip with an Internal Crack", Jut. J. 
Solids StAuotuAQ* , vol. 11, pp. 693-708, 1975. 


- 19 - 


14. S. G. Lekhnitskii, Theory of Elasticity of an Anisotropic Elastic 
Body , Holden-Day, Inc . , 196 kT. 

15. N. I. Huskhel ishvill , Singular Integral Equations , P. Noordhoff, 
Groningen, The Netherlands, 1953. 

16. F. Erdogan and G. D. Gupta, "On the Numerical Solution of Singular 
Integral Equations", £. App£. Moth., vol. 30, pp. 525-534 » 1972. 



r/A* 


APPENDIX A 


Expressions for the functions K-j and Kg and the constants (see 
equation 17) (material type I): 

cosh(s,ax) 

K l^ x,a ^ = Hjm^PTaT *- m 7 cosh(s^ah) ^ m l m l 0 B 5 ta nh(s 2 ah) 

cosh(s 2 ax) 

+ Vll e 5 ) - m a ^(3^ - ) - 

+ m^gpg tanhts^ah))] » 

-j cosh(s^ax) 

K 2 (x,a ^ = 2s ^/(aT ^ m 7 cosh(s^ah) ^ m 4 m ll P 5 

m-jn cosh(s 2 ax) 

" m 2 m 10 e 5 tanh ( s 2 ah ^ + m 8 cosh(s 2 ah) * 

m li 

♦ (jjT^ rn 2 ni 9^5 tanh(s-|ah) + nyi-,.,^)] , 

P(a) = ni 3 m 10 tanh(s 2 cih) - tanh(s 1 ah) , 

m, = l+v^s^j/Sg , m, = l« yx s 2 d z /6 5 , m 3 = 5,+^c, , 

m 4 “ s 2 +v yx c 2 ’ ra 5 = v xy +d lV S 5 ■ m 6 = v xy +d 2 s 2 /B 5 ’ 

”7 * Vl +C l ’ m 8 = V 2 +c 2 • m 9 = c l s r 1 

m 10 = C 2 S 2~^ y ^T) s ^l“ s l/^5 5 ^12 ~ c *2 _S 2^5 

m 13 = d r d 2 m ll /ni 12 * m 14 = ( m 5" m 6 I,1 n^ ,n 12^ 2m 13 ’ 



^15 = 1 ^15 K m 7^5^ m i m i o +m 4 in l ^ )/2s^ m^ 3 

m 17 - ■j"^n)^iHg)/2s*jHi'j 3 » 

n^| g c “niymi 1 pg (in^'i'H^nii g/ni-j 2 ) /2m-| 2^2 * 
m 19 = m 8 m l l b 5 Cm 3 +m g m 2 /m 7 2 )/2s 2 mi 3 . 


APPENDIX B 

Expressions for the function K 3 and the constants r.. (see equa 
tion 20) (material type II): 


K 3 (x,t,a) = — Q^ a j {[-r g sin(a) 2 ax)sinh(w 1 ax) 
1 A 


+ r 5 cos(w 2 ax)cosh(a)-|Cxx)] • [r 16 sin[ai 2 a(h~t)] 


• (r g sin(w 2 ah)cosh(w 1 ah) + r 1Q cos(to 2 ah) sinh(co 1 ah) ) 


- r 18 (^ cos[ l o 2 a(h-t)] - sin[w 2 a(h-t)])‘ 


•(r-j sin(w 2 ah)sinh(w 1 ah) + r 2 cos(w 2 ah) cosh^ah))] 


+ [r g sin(w 2 ax)sinh(w 1 ctx) + r fi cos(w 2 ax) cosh^ax)] 


• [-r 16 sin[> 2 a(h-t)] • (r g cos(w 2 ah)sinh(w 1 ah) 


- r 1Q sin(to 2 ah)cosh(i»i 1 ah)) + r 18 (^ cos0 2 a(h-t)] 



- sin[w 2 a(h-t)]) • (-r 2 sin(w 2 ah)sinh(iu^ah) 

+ cos(a> 2 ah) costi(w^ah))]} , 

Q(a) = r^g sinh(u^ah)cosh(w-|Cth) + r 2Q sin(w 2 ah)cos(u) 2 ah) . 
s-| = w-j+ i w 2 , s 2 = , uyo , 

C 1 = ^7 +1 * P 8 * c 2 = * d l = e 9 +1& 10 ’ d 2 = ^9 -1 ^l0 * 

r 1 = , r 2 = w 2 +v yx $ 8 , r 3 - l+Vy X (wj » 

r 4 “ v yx^ w 1^10 +w 2^9^ P 5 * r 5 = u l v xy + ^7 * r 6 = W 2 v xy +P 8 ’ 

r 7 = V xy + ^ W 1 P 9”“2 P 10^ P 5 * r 8 = *V v yx ’ r 9 = UJ l p 7" u 2 0 8' 1 * 
r 10 = W 2 & 7 +01 1% * r 11 = ”^10 +to 2 /f % ’ r 12 = p 9' u) 1^5 ’ 

r 13 “ "V B 10 r 12 /r ll » r 14 s "^ r 7 +r 8 r 12 /r ll^ 2r 13 » 

r 15 = r 4" r 3 r i2 /r n * r 16 = r 15 e 5 /4r 13 u l * r 17 ** “ r ir r lV r 11 

r 18 = Vl7 /4r 13 * r 19 ” r l r 10" r 2 r 9 ’ r 20 = r 2 r 10 +r i r 9 ’ 

r 21 * ” r 6 r 16 r 9" r 6 r 18 r l +r 5 r 10 r 16 +r 2 r 5 r 18 * 

r 22 ~ ^ r l r 6' r 2 r 5^ r 18V w 1 ’ 

r 23 = - r 6 r 10 r 16- r 2 r 6 r 18“ r 5 r 9 r 16- r l r 5 r 18 ’ 

r 24 = ^ r 2 r 6 +r l l "5^ r 1g (iJ 2 /w l * r 25 = _r 21“ r 24 * r 26 = r 21" r 24 

r 27 = r 22" r 23 ’ r 28 = r 22 +r 23 5 r 29 = “^l r 28 +w 2 r 26^^ w l +w 2^ 




